Practional-Power-Law Level-Statistics due to Dynamical Tunneling 
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For systems with a mixed phase space we demonstrate that dynamical tunneling universally leads 
to a fractional power law of the level-spacing distribution P{s) over a wide range of small spacings s. 
Going beyond Berry-Robnik statistics, we take into account that dynamical tunneling rates between 
the regular and the chaotic region vary over many orders of magnitude. This results in a prediction 
of P{s) which excellently describes the spectral data of the standard map. Moreover, we show that 
the power-law exponent is proportional to the effective Planck constant hag. 
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Spectra of quantum systems whose classical dynam- 
ics are either regular or chaotic usually show universal 
statistical properties. This fascinating relation between 
classical motion and quantum spectra is demonstrated in 
Ref . [l| , where it was argued that the spectra of generic 
regular systems show Poissonian statistics. In contrast, 
spectral correlations of classically chaotic systems can be 
described by random matrix theory [2,l3|. A justification 
of this conjecture was given in terms of periodic orbit 
theory [^-[6|. The nearest- neighbor level-spacing distri- 
bution P{s) is of central importance to the study of uni- 
versal spectral properties [3-Q . Results from these stud- 
ies are of broad interest for applications in, e.g. solid state 
physics [lOl, mesoscopic physics [III, cold atom physics 
[l^ . and atomic as well as acoustic physics |13j . 

The spacing distribution of generic Hamiltonian sys- 
tems has been the subject of an active debate over the 
last decades [14| - |22| . These systems show a mixed phase 
space, where disjoint regions of either regular or chaotic 
motion coexist [see Fig. [T]. Assuming statistically inde- 
pendent subspectra corresponding to regular and chaotic 
regions in phase space. Berry and Robnik computed the 
level-spacing distribution of mixed systems |14{ . In con- 
trast to the predicted level-clustering behavior, P{s) > 
at s = 0, numerically a fractional power-law distribution 

P{s) ex s^^ (1) 

for small spacings s with exponent /? G [0, 1] was observed 
(la [iq . Qualitatively this behavior may be described by 
the Brody distribution [2j, [2J] , as recently discussed in 
Refs. [1^ [ll|. Yet, this approach involves a free fitting 
parameter which is not related to any physical property 
of the system. 

Dynamical tunneling [171 |25I430| plays an important 
role for the level-spacing distribution, as it weakly cou- 
ples regular and chaotic states and thus enlarges small 
spacings between the corresponding levels. In Refs. |19l - 
|22| a phenomenological coupling strength between regu- 
lar and chaotic states was introduced, while in Refs. |17l . 
list a fit-free prediction of the level-spacing distribution 
was given. However, these results do not explain the nu- 
merically observed power-law distribution, Eq. ([1]). 

In this Letter we show that a fractional power-law dis- 



tribution universally arises over a wide range of small 
spacings because tunneling rates from different regular 
states range over many orders of magnitude. We give 
an analytical prediction of the level-spacing distribution, 
which is in excellent agreement with numerical data of a 
modified standard map [see Fig. [1] and a designed kicked 




FIG. 1. (color online) Level-spacing distribution P{s) at 
hcH = 1/100 of the standard map [3l|], with the classical 
phase space shown in the inset. We compare the numerical 
data [(green) histogram] to the analytical prediction, Eg. Q, 
(black solid line) and the Berry-Robnik result [1J| (black 
dashed line) (a) on a linear and (b) on a double-logarithmic 
scale. The typical couplings 2vm (triangles) mark the power- 
law regime [straight (red) line]. Below the smallest typical 
coupling 2iio one finds linear level repulsion. 



system [see Fig. [2]. Moreover, we demonstrate that the 
power-law exponent (3 scales like /3 oc /icff- For the small- 
est spacings below the power-law regime our result re- 
covers the well known linear level repulsion [8| . 

As model systems we study kicked systems described 
by the Hamiltonian H{q,p) = T{p) + V{q) J2nei. ^(* ~ "") 
[321 . By a stroboscopic view at integer times one gets 
an area preserving map on the two-torus. The rela- 
tive areas of the regular and chaotic regions in phase 
space are denoted by pr and pc = 1 — Pr- The map 
is quantum mechanically given by a unitary operator U 
on a Hilbcrt space of dimension N with effective Planck 
constant /icff = 1/iV. The semiclassical limit is ap- 
proached as ft-cff ^ 0. Solving the eigenvalue equation 
U\n) = e"^" \n) yields N eigenphases 0„ with eigenvec- 
tors \n) . According to the semiclassical eigenfunction 
hypothesis J33l - l35| one expects N,- « p^N regular states 
as well as N^ « p^N chaotic states. The unfolded spac- 
ings are s„ = (0n+i — 4'n)N/{2'K) for phases (f>„ which 
are ordered by increasing size. 

In order to study the influence of dynamical tunneling 
on spectral statistics we consider the standard map (36[ 
which is the paradigmatic kicked system. We use pa- 
rameter values for which it has one large regular island. 
To obtain significant statistics for small spacings we con- 
sider an ensemble of modified standard maps, where the 
regular island is preserved, while the chaotic dynamics is 
varied. Thus, the regular levels and the tunneling rates 
remain essentially unchanged, while the chaotic levels are 
strongly varied. This is realized by varying the kicking 
potential V{q) in the chaotic part [31|. Moreover, since 
partial barriers in the chaotic region may affect spectral 
statistics beyond dynamical tunneling [2Q|, we remove 
their influence by the above choice of the kicking poten- 
tial. Numerically, we focus on the /loff-regime with few 
regular levels (e.g. N,- = 6), as semiclassically {Ny -^ oo) 
the influence of tunneling on spectral statistics becomes 
less pronounced. We find a power-law distribution for 
P{s) at small s [see Fig. [1]. This behavior is also ob- 
served for another ensemble of kicked systems [331 [see 
Fig. [2]. 

We model the spectral statistics of systems with a 
mixed phase space by considering a random matrix 
Hamiltonian H which contains regular levels e^ and 
chaotic levels Sc on its diagonal. These levels are coupled 
by off-diagonal elements Vmj, which account for tunnel- 
ing contributions between the mth regular and the jth 
chaotic state. H is scaled such that the mean level spac- 
ing is unity and can be chosen real and symmetric for 
time-reversal invariant systems. The regular levels £r are 
semiclassically determined by the torus structure of the 
regular island [S^l- Since N^- is small, they do not show 
the Poissonian behavior assumed in Refs. [l^, [l^, [23 |. 
The chaotic levels Sc behave like eigenphases of a ran- 
dom matrix from the circular orthorgonal ensemble [7|. 
Here we assume that there are no additional phase-space 
structures within the chaotic region. The coupling ma- 
trix elements v^.j are modeled by independent Gaussian 



random variables with zero mean. The standard devia- 
tion Vrn of Vm,j docs not depend on the chaotic state j 
but is specific for each regular state m. In particular Vm 
is smaller for states m which quantize closer to the center 
of the regular island, vq < ... < w^Vr-i- The typical cou- 
pling Vm is related to the tunneling rate 7^ of the mth 
regular state by 



(2) 



which follows from the dimensionless form of Fermi's 
golden rule in kicked systems [331 • Hence, we model the 
probability density P{v) of all couplings by 
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(3) 



The tunneling rates 7^ are parameters of the random 
matrix model which can either be determined numeri- 
cally or analyticallY, e^. using the fictitious integrable 
system approach [28|, |30[ . Since the tunneling rates 7™ 
vary over many orders of magnitude, the typical cou- 
plings Vm embrace a wide range on a logarithmic scale 
[see the triangles in Fig. [T]. Hence, in contrast to previ- 
ous studies [il,[ia,[23 P{v) is not Gaussian but strongly 
peaked around small couplings. 

In the spirit of the semiclassical eigenfunction hypothe- 
sis [33I435I I we partition the level-spacing distribution into 
three distinct contributions 



P{s) =p,.„r(s) +Pc-c(s) +Pr-c(s). 



(4) 



Here, pr_i(s) describes the fraction of r-r spacings formed 
by two regular levels, Pc-c(s) the fraction of c-c spacings 
formed by two chaotic levels, and Pr_c(s) the fraction of 
r-c spacings formed by one regular and one chaotic level 
in the superposed spectrum [ij, [l^ [g^] ■ 

We evaluate the three contributions by making the fol- 
lowing assumptions: (i) The spacings of the chaotic sub- 
spectrum of H can be approximated by the Wigncr dis- 
tribution Pc{s) ~ 7rsp^/2c^'^'-*'''=' /^ with mean spacing 
1/Pc 0,8; (ii) consecutive regular levels are separated 
on scales larger than the mean level spacing. This is 
gencrically the case if there are less regular than chaotic 
states (pr < Pc) and hes is not much smaller than the reg- 
ular region pr (/leff ^ Pr)- Then one has just few regular 
levels, N-c w p^-fhes, which are semiclassically determined 
by the torus structure of the regular island. The inter- 
val between such consecutive regular levels then typically 
contains chaotic levels. 

The contribution of zeroth order r-c spacings s = {e^ — 
Ec I to the level-spacing distribution, neglecting couplings 
between regular and chaotic states, is given by 



Pr-c(s) == 2p,.pccxp 



'■n{spcf 



(5) 



Here Pc is the probability to have a chaotic level Eg in 
the distance s from the regular level Cr, J- Pc(s)ds = 



exp (— 7r(spc)^/4) is the probability to have no further 
chaotic level between e,- and Ec , and 2pi. is the probability 
of a zcroth order r-c spacing to contribute to P{s) [ij]. 

Dynamical tunneling leads to enlarged r-c spacings, 
which can be modeled by the 2x2 submatrices of H 



(6) 



This relies on degenerate perturbation theory [1^13 ^^.d 
is applicable because typically both r-r and c-c spacings 
are large compared to the couplings («„, ^ l/pr^c)- From 
Eq. © we calculate the tunneling improved r-c spacings 

s = a/s^ -I- 4z;2 such that 

oo oo 

p,_^{s) = fdv P[v) fdspi%{s) 6 (s - V52 + ivA . (7) 

—00 

This expression reflects that r-c spacings result from all 
possible zeroth order r-c spacings s and all possible cou- 
plings V, which are described by P{v). 

In order to compute the integrals in Eq. ([7]) 
we introduce polar coordinates {s'^ip) with s = 
s' cosif and 2v = s'sini^ [19|, such that p^-ds) = 



.(0) 



s J^ dtp P [^ sin if) pI J ^s cos ip). Calculating the re- 



maining integral gives 



Pr-c(s) 



(0) 



. Nr-1 ~ / 



with X (x) — y^n/2xe ^ I^Iq (a;^/4), where /q is the ze- 
roth order modified Bessel function of the first kind and 

Vm = Vral (l - 27rp2^^) . 

The contribution of r-r spacings to the level-spacing 
distribution is insignificant, due to assumption (ii) 

V.-M = 0. (9) 

The contribution of c-c spacings to P(s) is 



Pc-c(s) = Pc(s) [1 - PtS\ 



(10) 

where the first factor is the probability of finding a c-c 
spacing of size s in the chaotic subspectrum of H . The 
second factor describes the probability of having no reg- 
ular level within this c-c spacing, which is valid in the 
regime where s is smaller than all r-r spacings. 

Combining Eqs. ([8|)-(fT0| in Eq. (JH) gives our predic- 
tion of the level-spacing distribution in the presence of 
dynamical tunneling. Figures [1] and [2] show that this re- 
sult is in excellent agreement with spectral data of our 
example systems for s < 1. The small deviations for 
s > 1 can be attributed to approximation (i). 

Now we derive the fractional power- law of P(s), 
Eq. ([1]), in the tunneling regime. Here s is in between 
the smallest typical coupling ug and the largest typical 
coupling wat^-i. In this regime r-r spacings do not con- 
tribute [see Eq. ([9])]. Furthermore, the repulsion between 
chaotic levels [see Eq. ((TU)) ] ensures that c-c spacings only 
have significant probability for larger spacings and are in- 
significant deep in the tunneling regime (s ^ 1). Here, 







FIG. 2. (color online) Comparison of the level-spacing dis- 
tribution P(s) for the kicked system [331, with the classical 
phase space shown in the inset. The numerical data [(green) 
histograms], the analytical prediction, Eg. @, (black solid 
lines), and the Berry- Robnik result [IJj (black dashed line) 
are shown for h^a = 1/80 (lower curves), h^a = 1/120 (middle 
curves) and h^a = 1/160 (upper curves). The corresponding 
typical couplings 2vm (triangles, squares, and circles) mark 
the power-law regime [straight (red) lines]. The gray solid 
lines illustrate contributions from increasing v^ by increment- 
ing the upper summation index in Eq. Hill) from to A*'r — 1 
for /icff = 1/80. This gives rise to a staircase function. 



Eq. ([S]) reduces to p\_J^s) sa 2piPc which is used together 
! Wm in Eq. ([5]) leading to 



with u,, 



P{s) 



Af,-! 

N ^ 



2prP. 



X 



2v„ 



(11) 



X {s/2vj,i) behaves linearly for s < 2vjn and reaches a 
plateau of unit height for s > 2vm- Hence, one finds lin- 
ear level repulsion below the smallest coupling ug which 
would dominate the whole tunneling regime, if only one 
typical coupling was present in the system [iSl Il9l l22l | . 
However, the Um range over many orders of magnitude. 
This leads to an arrangement of the X functions, which 
individually behave like a step function, such that P{s) 
becomes a staircase on a double logarithmic scale [see 
Fig. [2]. This staircase resembles a power law, according 
to Eq. dD), and its slope is the power-law exponent 



/3 W log (iVr)/ log {vN,-l/vo). 



(12) 



This explicitly shows that the fractional power-law arises 
from typical couplings which range over many orders of 
magnitude. 

We now evaluate the scaling of the exponent /3 under 
variations of hcs. With Nr = Pr/heS, vn -1 ~ 1, and the 

I T I 1 

rough estimate vq ^ exp (— Cpr/^cff) [28|, l39|, we get 



/? (X /loff- 



(13) 



This type of scaling behavior, phenomenologically de- 
rived in Ref. [1^, is confirmed in Fig. [5] (hes = 1/80, 
1/120, and 1/160), with the almost constant ratio 
f3/hes = 26, 24, and 21 for the three cases of /leff- This 
result demonstrates the inapplicability of the Brody dis- 
tribution [2J,[23| for mixed systems, as it fails to simulta- 
neously describe the /loff-dependent power-law exponent 
for small spacings and the fixed Berry-Robnik-type dis- 
tribution of large spacings beyond the tunneling regime. 
Let us conclude by considering the level-spacing dis- 
tribution in the semiclassical limit h^g — )• 0. In this 
limit small r-r spacings appear such that assumption (ii) 
is violated. A generalized derivation shows that the r-c 
contribution to P{s) then still follows a power-law with 



exponent /? — > 0, as in Eq. ([T3|) . At the same time the 
r-r contribution approaches Poisson statistics. In combi- 
nation one thus recovers Berry-Robnik statistics in the 
semiclassical limit. 

To summarize, we have demonstrated how the wide 
range of dynamical tunneling rates universally leads to a 
power-law distribution of P{s) at small spacings s. We 
expect that this is also the fundamental mechanism which 
explains spacing-statistics in mixed systems with more 
complicated phase-space structures. 
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